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Since this is a meeting on inflation the title of this talk should probably have been "The cosmos as a 
lab for inflation" . I will illustrate how, by combining WMAP observations of the cosmic microwave 
background (CMB) with smaller scales CMB experiments and large-scale structure surveys, it is 
possible to constrain inflationary models. WMAP data in combination with these external data 
sets, offers a window to the very early universe and the physics of inflation: enables us to confirm 
the basic tenets of the inflationary paradigm and begin to quantitatively test inflationary models. 
Understanding the complex physics and complicated systematics of the low redshift observations is 
crucial to make progress in this direction. 



I. INTRODUCTION 



The first year results from the Wilkinson Microwave 
Anisotropy probe (WMAP) were announced on Febru- 
ary 11 2003. On the same day the telescope was re- 
named in honor of Prof. David Wilkinson, member of 
the science team and pioneer in the study of cosmic mi- 
crowave background (CMB) radiation. The primary goal 
of the WMAP mission is to produce a high-fidelity all- 
sky polarization-sensitive map of the cosmic microwave 
background (CMB) radiation to determine the cosmol- 
ogy of our Universe. After a year of observation it has 
produced a full sky map of the microwave sky in 5 fre- 
quencies, with a resolution a factor 30 higher than the 
previous full sky map as produced by the COBE satellite 
0. This is the cleanest picture of the early Universe; 
the structures on the CMB -the pattern of hot and cold 
spots- carry information about the composition, geom- 
etry, age etc. of our Universe (e.g., 0, El). Here I 
will illustrate how it is possible to constrain inflation- 
ary models by combining WMAP's CMB observations 
with external data sets: in particular, the smaller scales 
CMB experiments CBI and ACBAR, and observations of 
large-scale structure probes such as the power spectrum 
of Two-degree-field galaxy redshift survey (2dFGRS; 0) 
and the linear matter power spectrum as recovered by 
PI and from Lya forest observations. We also con- 
sider the Type 1A supernoave results |26l_and HST Key 
project Hubble constant determinations M . The results 
of this analysis were reported in [25l l28l l30| . 



*The WMAP science team composed by: C. Barnes, C. Bennett 
(PI). M. Halpern, R. Hill, G. Hinshaw, N. Jarosik, A. Kogut, E. 
Komatsu, M. Limon, S. Meyer, N. Odegard, L. Page, H. Peiris, D. 
Spergel, G. Tucker, L. Verde, J. Weiland, E. Wollack, E. Wright. 



II. WMAP 

If the CMB is Gaussian (see [13) then the informa- 
tion in the mega-pixel maps can be compressed in the 
power spectrum [I^. We show the WMAP tempera- 
ture power spectrum in Fig. 1 (left), where the points 
with error-bars are the (band power) data with the error 
given by the instrumental noise. The solid line is our 
best fit model and the gray area shows the cosmic vari- 
ance uncertainty. At the moment we can reliably measure 
only the temperature-polarization (TE) cross-correlation 
power spectrum which is shown in Fig. 1 (right) where 
the points with error-bars are the band-power data. The 
temperature power spectrum makes precise predictions 
for the TE power spectrum on large scales. The solid 
line in Fig. 1 (right) is this prediction from the best fit to 
the TT data only. For £ > 10 the data are in remarkable 
agreement with the TT predictions, this is a triumph for 
the standard cosmological model (see H^); the WMAP 
detection of TE anti-correlation at I ~ 50 — 150 rules out 
a broad class of active models and implies the existence of 
super-horizon adiabatic fluctuations at decoupling: both 
inflation and the Ekpyrotic scenario predict the existence 
of super horizon fluctuations. 

For I < 10 we detect a signal that cannot be seen in 
the TT data alone: this is the signal of early re-ionization 
Such a strong signal at low multipoles gives a con- 
straint on the integrated optical depth to the last scatter- 
ing surface r = 0.17±0.04, for instantaneous reionization 
this corresponds to a reionization redshift of z re = 17 ±3. 

Why WMAP has such a good signal to noise? A part 
from the fact that it is an all sky experiment which ob- 
served the sky for a full year, this is the result of an 
obsession of all those involved in designing and building 
the satellite for minimizing systematics and maximizing 
the signal to noise. During the designing and building 
of the satellite the emphasis was on minimizing system- 
atic errors so that the data could be straightforwardly 
interpreted. For example there are 10 different channels 
from 20 to 95 GHz, so that foreground subtraction is not 
a problem [il. Il2|. Undesirable 1// noise is minimized 
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by the design of WMAP radiometers 0], uncorrelated 
noise between pixels, accurate determination of the beam 
patterns 0, 122] and well understood properties of the ra- 
diometers are invaluable for this analysis. 

The rest of the analysis must therefore have the same 
level of "obsession" . For example, to extract as much in- 
formation as possible from the maps without introducing 
systematics, the TT power spectrum is obtained by op- 
timally combining "only" 28 cross-channels correlations. 
This method has then be checked against another 2 in- 
dependent methods [l^. Also, the point source subtrac- 
tion has been done in 3 different ways HE 03 (they all 
agree!). Finally in the covariance matrix we propagate 
all the effects we know about (beam errors, noise, sky 
cut etc.) To cut a long story short, we spent a long time 
worrying about 2% error on the error! There are effects 
that we are aware of, which we neglect; these effects give 
about 0.5 to 1% error on the error. 
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FIG. 1: Left: The TT power spectrum obtained from an op- 
timal combination of 28 cross-power spectra. The solid line is 
the best fit model, and the gray area shows the cosmic vari- 
ance error. Right: the temperature-polarization cross corre- 
lation power spectrum. The solid line is the prediction for the 
polarization signal from the TT data. The excess power at 
£ < 10 is due to early re-ionization (see |T^|h 



sian, then the spherical harmonics coefficients ai m are 
normally distributed, but the angular power spectrum, 
Ci, is not. For large I the Central Limit Theorem will 
make the Gaussian likelihood an increasingly good ap- 
proximation, but the error-bars in the WMAP data are 
so small, that previous likelihood approximations exten- 
sively used in the literature are not suitable. We intro- 
duce a different form of the likelihood function that is 
accurate to better than 0.1%. The likelihood function 
and its covariance matrix has been calibrated on 100,000 
Monte Carlo realization of the sky as seen by WMAP 
|30|. This calibration is also important because we can 
now compute the reduced % 2 for the "best fit model" and 
use it to compute the probability that the temperature 
pattern observed in the sky is a realization of the best fit 
model. 



We perform a likelihood analysis for the cosmological 
parameters: at each point in parameter space (i.e., for 
a given set of cosmological parameters) a new model Ce 
is computed using the publicly available code CMBfast 
4.1 [23|. Even with a supercomputer, at 2 seconds per 
evaluation, with a naive grid-based approach it would 
take a prohibitive amount of CPU time to perform the 
analysis for 6 or more parameters. We resort to Markov 
Chain Monte Carlo (MCMC): this approach has become 
the standard tool for CMB analyses jEHEiil. In fact, 
for a flat reionized universe we can evaluate the likelihood 
~ 120, 000 times in < 2 days: this is adequate for finding 
the best fit model and reconstruction of the 1- and 2- a 
confidence levels for the cosmological parameters. While 
more details are in |30| and references therein, here we 
just remind the reader that a MCMC is similar to a ran- 
dom walk in parameter space: at each step a new point 
is sampled in parameter space but the density of steps 
(points in parameter space) is proportional to the a pos- 
teriori likelihood. After an initial "burn-in" period, the 
chain has no dependence on the starting location. The 
important point to keep in mind about MCMC is that a 
MCMC cannot tell us about regions of parameter space 
where it has not been: it is very important that the chain 
achieves good mixing (i.e., moves rapidly in parameter 
space) and fully explore the likelihood surface. Using an 
MCMC that has not fully explored the likelihood sur- 
face for determining cosmological parameters will yield 
wrong results, thus it is important to have a convergence 
criterion and mixing diagnostic. For any analysis of the 
WMAP data we strongly encourage the use of a conver- 
gence criterion, for example we use the method proposed 
by M- 



III. THE METHOD 

The analysis method must match the level of accu- 
racy and precision of the data. If the CMB sky is Gaus- 



We speed up convergence and improve mixing by using 
the physical parameters of 3 and [l8| (e.g., physical cold 
dark matter density Q, c h 2 instead of Q, c and the charac- 
teristic angular scale of the acoustic peaks @a instead of 
the Hubble constant). 
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A. WMAPonly results 

We first present the cosmological parameter con- 
straints from WMAP data only. Then we consider exter- 
nal data sets (smaller scale CMB experiments and large 
scale structure power spectrum measurements) and com- 
pare these data to the predictions of the best fit model for 
WMAP data. We find that this extrapolation of WMAP 
predictions to low redshift and smaller scales works re- 
markably well. We thus proceed to perform a joint like- 
lihood analysis of WMAP data with these external data 
sets. Apart from the fact that the CMB seems to be 
Gaussian |17| . consistent with the generic prediction of 
inflationary models, and that the Universe seems to have 
been reionized much earlier that previously though |16| . 
the most striking result, probably, is that the "vanilla", 
flat LCDM model still fits the data amazingly well: 6 pa- 
rameters fit 1348 data-points (these are the power spec- 
trum multipoles, one could be bolder and say that 6 pa- 
rameters fit ~ 10 6 data, these being the pixels!). The 
parameters that best fit WMAP data are summarized in 
table 1. 

Remarkably, this same model fits not only WMAP 
data but also other small scale CMB experiments and a 
host of other cosmological observations: the Hubble con- 
stant constraint (Hn = 72 ± 5) is in agreement with the 
Hubble key project [!j determination of H = 72 ± 3 ± 7 
at the 68% confidence level(C. L.); the age determination 
of 13.4 ± O.SGyr is in good agreement with other age- 
determination based on stellar ages; the baryon abun- 
dance flbh = 0.024 ± 0.001 implies a primordial deu- 
terium abundance in good agreement with D/H determi- 
nation from quasar absorption lines; the determination of 
amplitude of fluctuations parameterized by as (the r.m.s. 
fluctuations smoothed on a scale of 8Mpc/h) is in good 
agreement with clusters and weak lensing estimates. 



IV. EXTERNAL DATA SETS 

We also consider smaller scales CMB experiments 
ACBAR [ig| and CBI [Hill on scales where they do not 
overlap with WMAP measurements and they are not af- 
fected by possible SZ signal. We combine CMB data with 
measurement of the low redshift universe. Galaxy red- 
shift surveys allow us to measure the galaxy power spec- 
trum at z ~ and observations of the Lyman a (Lya) 
forest allows us to probe the dark matter power spec- 
trum at z ~ 3. We use the Anglo- Australian Telescope 
Two Degree Field Galaxy redshift Survey (2dFGRS;0), 
and the linear matter power spectrum as recovered by 
0| and from Lya forest observations. When analyz- 
ing the 2dFGRS survey we find that we need to model 
the non-linear evolution of the power spectrum and non- 
linear redshift space distortions, even when limiting the 
analysis to a regime traditionally considered to be still 
well described by linear theory. We tested our modeling 
extensively on mock catalogs of the 2dFGRS. The goal 



of our modeling is to relate the shape and the amplitude 
of the observed power spectrum to that of the under- 
ling linear dark matter power spectrum as constrained 
by CMB data. Thus we use the information about the 
redshift space-distortion parameters and galaxy bias as 
in [24[ and [2|j (for more details see |30] ) . 

For the matter power spectrum derived from Lya for- 
est observations we follow the approach on [Tl| and in- 
crease the error-bars where their analysis is not in agree- 
ment with that of ||). We stress here that the addition 
of Lya forest data appear to confirm trends seen in the 
other data sets and tighten cosmological constraints how- 
ever our results do not rely on Lya data. More observa- 
tional and theoretical work is still needed to confirm the 
emerging consensus that the Lya forest data traces the 
large scale structure (see L. Hui and U. Seljack talks) 
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FIG. 2: The data set, and the LCDM model best fit to 
WMAP data. In the top panel WMAP data (black) CBI (red) 
ACBAR (blue) where the multipole £ has been converted to 
wavenumber k. In the bottom panel 2dFGRS (black) and Lya 
(green). The extrapolation of this model to low redshift and 
to smaller scales fits the data remarkably well. For 2dFGRS 
we have used redshift space distortions parameters: /3 = 0.45, 
a v = 380 km/s, and bias parameter b = 1.13. 

The key cosmological parameters that are relevant 
for single field inflation arc: the constraint on the to- 
tal matter-energy density of the universe relative to the 
critical density il tot , the spectral slope of the primor- 
dial power spectrum n s the running of the spectral index 
dn/dlnk and the tensor to scalar ratio r. These are re- 
ported in table 1 and in figure . 

V. RESULTS AND CONCLUSIONS 

I will present here 2 sets of conclusions: 

• 1) for the pessimists I have shown how to pollute 
the clean CMB data with the messy gastrophysics 
of the low redshift probes of large-scale structure. 
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FIG. 3: The parameter space spanned by viable slow-roll inflation models, with the WMAP+CBI+ACBAR+2dFGRS+Lya 
68% confidence region in dark blue and 95% confidence region in light blue (see (25| for more deatils). 



• 2) for the optimists, I have shown that all these 
different and independent data sets are in remark- 
able agreement. Within the context of the inflation 
+ LCDM paradigm cosmological parameters con- 
straints obtained from these different datasets are 
consistent. The results obtained at z ~ 1088 can 
be extrapolated forward to z ~ and describe a 
host of astronomical observations. 

This gives us the confidence to extrapolate our re- 
sults backwards, to well before z ~ 1088 and place 
constraints on Inflation, as presented in j2f| (see 
also D. Spergel talk and fig. 3). 

For both optimists and pessimists I hope I have illus- 
trated how it is possible to enhance the scientific value 
of CMB data from z ~ 1088 by combining it with mea- 
surements of the low redshift universe. When including 
these external data sets one should keep in mind that 
the underlying physics for these data sets is much more 
complicated and less well understood than for CMB data 
and especially WMAP data. For both pessimists and 
optimists I hope I have highlighted the importance of 
improving the rigor of large-scale structure data analy- 
sis and of the interpretation of these observations, and 
how crucial it is to quantify an d account for the level of 
systematic contamination. 

The data set, the WMAP temperature and temperature- 
polarization angular power spectra are available at 



://lambda.gsfc. nasa.gov If using these data set 
please refer to Hinshaw et al. (2003) for the TT power spec- 
trum and Kogut et al. (2003) for the TE power spectrum. 
We also make available a subroutine that computes the 
likelihood for WMAP data given a set of Ce, please refer to 
Verde et al. (2003) if using the routine. 
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